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Figure 2. Near-infrared contour emission map obtained for the H2 1-0 S(l) transition on the HH46/47 object. The position of 
the HH objects are labelled in this figure, as well as the bright and extended counter-jet flow. The spatial co-ordinate is given 
in IRSPEC rows (each corresponding to a length of 2.2" per pixel on the sky) on the left axis. To plot this image, we used 
7 contours increasing in steps of \/2 from a base level of 0.25 X 10 -11 ergs _1 cm~ 2 (im _1 . The vertical line is plotted at zero 
heliocentric radial velocity. 
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Figure 3. Position-velocity emission maps obtained for the higher excitation H2 line emission discussed in the main 
text. The 2-lS(l) transition has contours plotted linearly from a base contour of 0.09 X 10 — 1:L ergs _1 cm~ 2 /im — 1 , and 
stepped by 0.05 X 10 _11 crgs _1 cm~ 2 (im _1 . The 1-0S(0) transition has contours plotted linearly from a base contour of 
0.13 X 10 — 11 ergs _1 cm~ 2 /im -1 , and stepped by 0.07 X 10~ 11 ergs _1 cm - 2 (im _1 . The 2-lS(2) transition has contours plot- 
ted linearly from a base contour of 0.11 X 10 _11 crgs _1 cm~ 2 /im -1 , and stepped by 0.025 X 10 _11 ergs -1 cm~ 2 /im -1 . The 
3-2 S(3) transition has contours plotted linearly from a base contour of 0.16 X 10 _11 erg s _1 cm~ 2 /im -1 , and stepped by 
0.02 X 10~ 11 ergs~ 1 cm -2 /jm" 1 . 
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Figure 4. H2 excitation diagrams for the counter-jet. The derredened column densities are plotted against the transition 
upper level energy for the four H2 lines (labelled next to the plotted data point) observed with the highest signal-to-noise (see 
Table 3). The arrow in the plot represents a 3<r upper limit. The solid line represents the least squares fit to the errorless data 
points providing an excitation temperature of T ex = 2600 K (± 300 K) for the whole counter-jet (Total; bottom diagram) and 
the northern peak in the counter-jet (North; top diagram) 
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Figure 5. Position- velocity emission maps obtained in the [Fell] 1.257 /im and [Fell] 1.644 fim transitions. For the [Fell] 
1.644 (im, the base contour is plotted at 0.55 X 10 _11 ergs _1 cm~ 2 /im _1 with a linear step of 0.15 X 10 —11 ergs — 1 cm -2 /jm" 1 
whereas for the [Fen] 1.257 (im, the base contour is at 1.1 X 10 _11 crgs _1 cm~ 2 ^m _1 and adjacent contours are stepped 
linearly by 0.7 X 10 _11 crgs _1 cm~ 2 (im -1 . In both maps, the vertical solid lines are drawn at zero velocity emission whereas 
the horizontal solid lines is drawn at the position of the continuum emission from the HH46 object. These positions have been 
calculated using the spatial offsets given in Table 1. 
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Figure 7. Overlay of the H2 1-0 S(l) and [Fen] 1.644 /im emission maps. 
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ABSTRACT 

Long slit spectra of the molecular outflow Herbig-Haro (HH) 46/47 has been 
taken in the J and K near-infrared bands. The observed H2 line emission confirms the 
existence of a bright and extended redshifted counter-jet outflow southwest of HH46. 
We show that this outflow seems to be composed of two different emission regions 
which have distinct heliocentric velocities, in contrast with the optical appearance, 
and which implies an acceleration of the counter-jet. 

The observed [Fe 11] emission suggests an average extinction of 7-9 visual magni- 
tudes for the region associated with the counter-jet. 

Through position- velocity diagrams, wc show the existence of different morpholo- 
gies for the H2 and [Fe 11] emission regions in the northern part of the HH 46/47 outflow. 
We have detected for the first time high velocity (-250 km s -1 ) [Fen] emission in the 
region bridging HH46 to HH47A. The two strong peaks detected can be identified 
with the optical positions B8 and HH47B. 

The H2 excitation diagrams for the counter-jet shock suggest an excitation tem- 
perature for the gas of T ox 2600. The lack of emission from the higher energy H2 
lines such as the 4-3 S(3) transition, suggests a thermal excitation scenario for the 
origin of the observed emission. Comparison of the H2 line ratios to various shock 
models yielded useful constraints about the geometry and type of these shocks. Pla- 
nar shocks can be ruled out whereas curved or bow-shocks (both J and C-type) can 
be parametrised to adjust our data. 

Key words: ISM: individual objects: HH46/47 - ISM: jets and outflows - ISM: 
kinematics and dynamics - Infrared: ISM: lines and bands 



1 INTRODUCTION 

For its tight spatial alignment of the emission knots and its 
clear visibility against the progenitor dark globule, HH 46/47 
is one of the most remarkable examples of the classical 
Herbig-Haro object. Being optically bright, this HH com- 
plex has been the subject of many studies of its spectral, 
photometric and kinematical properties after its discovery 
by Schwartz (1977). 

The morphology of the Gum nebula region (Dopita 
et al. 1982) shows the existence of a very collimated jet that 
is composed of several components which extend about 0.39 
pc on the sky at the assumed distance of 350 pc (Eisloffel 
and Mundt 1994). A young T Tauri star still obscured by 

* Based on observations collected at the European Southern Ob- 
servatory, La Silla, Chile 



the surrounding dust envelope sits at the center of this sys- 
tem. This star illuminates a cloud of dust revealing a re- 
flection nebula from which the optically visible HH46 and 
HH 47(A,B,C,D) jet components emerge. 

The bipolar nature of this outflow was first ascertained 
by Graham and Elias (1983), on discovering blue and red- 
shifted outflowing shocked gas to either side of the central 
source. Schwartz et al. (1984) reported the first proper mo- 
tions for some dense knots along the flow suggesting the 
occurrence of episodic ejections of gas material from the en- 
ergy source. Combining the radial velocity measurements 
(Meaburn & Dyson 1987, Hartigan et al. 1993, Morse et al 
1994) with a detailed study of the proper motions of the 
knots in the HH 46/47 jet (Eisloffel & Mundt 1994), a space 
velocity of about 300kms _1 is found for the jet. 

While several imaging, spectroscopic and proper mo- 
tion studies have been carried out in the optical for the 
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blueshifted part of HH 46/47, there has been considerably 
less work done for the redshifted part of the flow. This 
well-collimated counter-jet flow located southwest of the 
source was discovered by Reipurth and Heathcote (1991) in 
[SII] AA6717/6731 line emission and found to be redshifted 
with a heliocentric radial velocity of about +125kms _1 . PC- 
band imaging of the redshifted counter-flow has been per- 
formed by Eisloffel et al. (1994), together with H 2 1-0 S(l) 
(2.1218 Aim) and [Fell] (1.644 Aim) line imaging. More re- 
cently, proper motion studies of the H2 emission have been 
carried out for the H2 knots in HH 46/47 (Micono et al. 
1998). 



2 OBSERVATIONS AND DATA REDUCTION 

The HH 46/47 complex was observed in 1995 March 14 to 
17 using the cryogenically cooled grating spectrometer IR- 
SPEC installed on the ESO 3.5 m New Technology Telescope 
(NTT) sited at La Silla, Chile. The instrument was equipped 
with an InSb CCD detector comprising 58 x 62 square pix- 
els which provided a pixel scale of 2.2 arcseconds on the 
sky. We have positioned the 2 arcminute long x 2.2 arc- 
seconds wide slit at a position angle of 57 degrees over the 
HH 46/47 bipolar flow, centered on approximately the HH 46 
emission nebula continuum (Figure 1). This orientation on 
the sky coincides with the outflow axis of both optical and 
infrared maps of HH 46/47 and the 2 arcminute long slit 
allowed us to measure simultaneously the emission region 
from the HH47A object to the counter-jet flow. The precise 
positioning of the observed slit position was set by peaking 
up the H2 emission maximum. The telescope pointing accu- 
racy, normally better than 1" (Moorwood et al. 1991) has 
been checked with observations of the standard stars. 

The seeing disk remained under 1" with an average 
FWHM of 0.8". Grating no. 2 was used with the slit set to 
a width of 2 pixels. The spectral resolution was ~ 0.001 fim 
(FWHM) at 2.122 fim but varied little over the whole K 
band region (R — A/ A A ~ 2300). For the two observed 
[Fell] transitions, the resolution was R « 3300 and 1600 at 
1.257 /im and 1.644 Aim, respectively. At these resolutions, 
every frame contained a very short wavelength range and 
only one emission line was registered per frame. Because 
of the extent of the HH 46/47 outflow on the 2' slit, ex- 
plicit sky observations were taken by offsetting the tele- 
scope 40" south from the nominal object position. The 
integration times were 60 seconds on source and 60 sec- 
onds on sky and a final frame was determined by averaging 
many integrations. The observational procedure consisted in 
taking several exposures at the different programme wave- 
lengths on each night and then coadding them to produce 
a higher signal-to-noise final frame. The observations were 
performed on the following emission lines of H2 arising in 
the K band; v=l-0S(l) [2.1218 ^m], 1-0 S(0) [2.2235 Aim], 2- 
1S(1) [2.2477 Aim], 2-1 S (2) [2.1540 fim] , 2-1 S(3) [2.0650 Aim], 
3-2 S(l) [2.3846 Aim], 3-2 S(3) [2.2014 Aim], and the 4-3 S(3) 
[2.3445 Aim]. In addition, we have observed the Br7 tran- 
sition at 2.166 Aim and two emission lines of iron: [Fell] 
1.257 Aim and [Fen] 1.644 /im arising on the J and H bands, 
respectively. Table 1 provides a list of all the line observa- 
tions together with the observational parameters for each 
line transition. 



One of the main goals of the project was to observe 
the high-excitation H2 lines (with upper level energies above 
13000 K) of molecular hydrogen in the K band, especially 
the ones arising from vibrational levels v'=3 (e.g. the 3- 
2S(3) [2.2014 A«m]) and v'=4 (e.g. the 4-3 S(3) [2.3445 Aim]). 
Making exposures at such low-light levels that are character- 
istic of these faint lines in shocked regions required a special 
technique since several effects such as the non-linearity of the 
array at low-light levels may be comparable with the photon 
noise level (Gredel and Moorwood 1991). The method em- 
ployed consisted in taking several frames at the same wave- 
length but slightly shifted by 1/2 and 1 pixel, thus ensuring 
that the same flux is sampled by different pixels on the ar- 
ray. This also helps to prevent the spread of bad pixels to 
the final image. To produce the final reduced image we align 
the frames to a common wavelength grid and coadd up the 
exposures to gain sufficient signal-to-noise on the emission 
line. 

Three photometric standard stars were routinely ob- 
served throughout the night at the same object wavelengths 
and airmass to allow a proper cancellation of telluric absorp- 
tion features and to flux calibrate the object. The stars ob- 
served were: Beta Car (HR 3685), Eta Cma (HR 2827) and 
Iota Car (HR 3699). The relative variation of the standard 
star flux over the three night period, provided an estimate 
of the uncertainty on the calculated flux of ~ 10%. After 
inspecting the reduced standard star frames, we noted that 
at two wavelengths, the standard star continuum displayed 
deep absorption features at 2.166 Aim and 1.642 Atm. These 
features are known to be associated with the absorption by 
the stellar photosphere of the hydrogen recombination ra- 
diation. The lines in question are the Hi [n=7-4] and Hi 
[n=12-4], for the 2.166 Aim and 1.642 Aim, respectively and 
were artificially removed by linear interpolation between the 
line edges, a common practice since the error committed is 
negligible compared to the star's continuum absolute flux. 

The flux calibration procedure required a flux table for 
each standard containing a properly sampled representation 
(lambda vs. flux) of the continuum for each standard star. 
Such file was constructed for each star based on the broad 
band flux given in the SIMBADp] database. Using this table, 
it is possible to reproduce the star's flux distribution over 
the wavelength region of interest and to use this curve to 
calibrate the flux counts in the object frames. All data were 
primarily reduced with specific ESO routines available in 
MIDAS for IRSPEC observations. The later stages of data re- 
duction that include standard procedures (sky subtraction, 
flat-fielding, wavelength and flux calibration) were assisted 
with routines from the STARLINK project. 

With the current IRSPEC standard setup, a slight spa- 
tial displacement of the objects occurs between frames taken 
at different wavelength settings. This is known to be caused 
by changing the grating settings between the different expo- 
sures (Gredel and Moorwood 1991). Thus, in order to correct 
for this effect, we have determined the relative spatial shift 
between the different wavelengths frames by calculating the 
centroid position in the spatial direction for the observed 
star continuum. This was achieved by fitting a gaussian (in- 

I The SIMBAD database is operated by the Centre de Donnees 
astronomiques de Strasbourg, France 
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strumental spatial profile) to the spatial profile of the star's 
image for all wavelengths. We found a shift of +3.5 pixels 
for the [Fell] 1.257pm image map and of —2.5 pixels for the 
[Fe n] 1.644 /im image map relatively to the Hb 1-0 S(l) map. 
All the other H2 lines have intermediate shifts ranging from 
-0.5 at 2.065 pm to +2.5 at 2.3846 pm. All the calculated 
spatial shifts between the frames are shown in column 5 of 
Table 1. 



3 RESULTS 

A total of eight H2 transitions were observed during the first 
two nights. These arise from vibrational levels: v=l, 2, 3 and 
4, and have all rotational jumps of A J = +2. Table 1 shows 
the parameter list of the observed emission lines. For each 
line observed, this table shows in the fourth column the total 
exposure time of the final reduced frame. In addition, the 
fifth and sixth columns give respectively, the spatial shift 
used to correct the different grating settings to the same 
offset and the velocity resolution for each frame. 

The heliocentric velocity correction for the time the ob- 
servations were performed (— 3kms _1 ) is negligible when 
compared with the width of a pixel (48-98 km s _1 ). Thus, 
we have chosen not to make any corrections on the calcu- 
lated position- velocity maps after converting the wavelength 
scale to the radial velocity scale on the maps. This conver- 
sion was applied only to the H2 1-0 S(l) and the two [Fell] 
maps since only these are bright enough to compute the 
radial velocity of emitting gas. The velocities are thus mea- 
sured with respect to the Bok globule which is believed to 
have a heliocentric radial velocity ~ +22kms~ 1 (Hartigan 
et al. 1993). 



3.1 The line emission 

Figure 2 shows the position-velocity contour emission map 
observed for the H 2 1-0 S(l) transition at 2.1218 pm. This 
transition is one of the strongest emission lines of H2 to 
arise in the K band and has been extensively used to map 
regions of shocked gas in molecular clouds (see for example 
Eisloffel 1997, for a recent review). Some unwanted features 
are also present in this emission map; some CCD columns 
to the right and left of the main central object chain show 
vertical and narrow linear emission structures that stand 
above the background noise level. These are probably due 
to the combined effect of poor cancellation of sky lines (due 
to the intrinsic sky emission line intensity changing over 
the integration period) in the reduction procedure and to 
the slight movements of the grating after changing (wave- 
length) position. The wavelength position of these features 
were compared with the list of OH airglow emission lines 
given in Olivia and Origlia (1992) and can be positively 
identified with the residual sky emission lines left over from 
the reduction procedure. This allows to check the on-line 
(mechanical) wavelength calibration (usually precise to ~ 1 
pixel). Using the two brightest OH lines in our spectra at 
2.11710pm [(9,7) Ri(2)] and 2.12440pm [(9,7) fli(l)], we 
have calculated the deviation of their profile centroids to 
the ones listed in Oliva & Origlia (1992). This gives a r.m.s. 
deviation error of 0.46 pixels and therefore imply an asso- 



ciated uncertainty of « 30kms _1 for the measured radial 
velocities. 

Apart from the continuum emission that is associated 
with the HH46 reflection nebula near the infrared source 
located close to the middle of the frame (rows ~22-25), 
there are several objects clearly seen as bright emission knots 
aligned vertically in the middle of the array. These objects 
are labeled in Figure 2 and their relative positions can be 
identified using the maps of Eisloffel et al. (1994) and Mi- 
cono et al. (1998). 

Figure 2 shows the northeastern blueshifted outflow to 
be mainly composed of two bright emission knots, HH46 
and HH47A. In the middle of the frame, HH46 appears as 
a bright knot blended with the nebular continuum emis- 
sion from the "reflection nebula" centered at ~ row 24. The 
emission from HH46 appears blueshifted with a heliocentric 
radial velocity of about — 70kms _1 , which is at the lower 
end of the radial velocities measured for [SII] and Ha emis- 
sion (Dopita et al. 1982, Meaburn & Dyson 1987, Reipurth 
& Heathcote 1991, Hartigan et al. 1993, Morse et al. 1994) 
which lies in the range w —80 to — 170 km s -1 . 

Further north, at the top of the PV map, there is a 
bright knot which we associate with the infrared counter- 
part of the high-excitation HH47A object. This object is 
commonly associated with the brightest working surface of 
the northeastern outflow and is located at a distance of 
about 65" northward of HH46. The centroid velocity of this 
bright emission peak is seen at a heliocentric radial veloc- 
ity of about — 100 km s -1 , which correlates well with the 
position of the optical, Ha and [SII], emission knots seen 
blueshifted at ~ -80 to -140 kms -1 (Dopita et al. 1982, 
Meaburn & Dyson 1987, Hartigan et al. 1993, Morse et 
al. 1994). HH47A is considered the best example of a bow 
shock/Mach disk typical interaction, even though it is a sec- 
ondary ejection from the source that is ploughing through a 
pre-accelerated and dense medium - the wake of an earlier 
shock now observed as HH47D. Imaging the HH47A shock 
in the near-infrared H2 1-0 S(l) line shows that it splits up 
in two regions: a bright knot Z and a fainter nebulosity to 
the south named as knot Y (Micono et al. 1998). In the 
region bridging HH46 to HH47A, there appears to be no 
significant H2 emission in our maps, apart from the faint 
emission filament located north of HH46 and extending up 
to a distance of 20" (rows ~ 25-32) at essentialy the same 
radial velocity. This is in good agreement with near-infrared 
maps; indeed, only a fainter emission patch can be seen to 
the west of the knot Z (Micono et al. 1998), which lies out- 
side our slit. Even though this object lies at the edge of the 
frame, where the detector array may suffer from vignetting, 
we have, nevertheless, measured the line intensity and es- 
timated the integrated column density for the brightest H2 
lines in this object (see Table 3) . 

The most striking feature of the H2 near-infrared maps, 
however, is the bright and extended emission at the base 
of the redshifted outflow located south of HH46. This emis- 
sion is associated with an extensive counter-jet shock clearly 
visible to the southwest of the HH 46/47 source in the H2 
1-0 S(l) image presented by Eisloffel et al. (1994). This 
highly collimated counter-jet is not detected in Ha emis- 
sion and was only discovered by Reipurth and Heathcote 
(1991) in their [SII] AA6716/6731A emission map as being 
composed of three emission knots extending 11" on the sky. 
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The counter-jet faint appearance in the optical is probably 
due to the increased extinction at optical wavelengths and 
the fact that the counter-flow is moving towards the interior 
of the Bok globule. However, it is readily seen as the bright- 
est object in our near-infrared PV maps since the extinction 
decreases about an order of magnitude from optical to in- 
frared wavelenghts. This counter-jet is known to be at the 
base of a larger (infrared) bow-shock that extends deeper 
into the interior of the Bok globule for about 2' until it es- 
capes out of the globule at HH47C (also seen in the infrared; 
e.g. Eisloffel et al. 1994). Figure 2 shows the counter-jet to 
be resolved in to two spatially blended emission knots but 
that appear clearly displaced in heliocentric radial velocity. 
The extent of this counter-jet is 22" on the sky or 0.037 pc 
at the distance of HH 46/47 (using a distance to the globule 
of 350 pc, Micono et al. 1998). The peak-to-peak separation 
of these knots is 0.015 pc (8.8" on the sky) which correcting 
for projection effects of the outflow orientation angle with 
respect to the plane of the sky (~35 degrees, Eisloffel and 
Mundt 1994) gives a separation distance of 0.026 pc. 

The most interesting feature seen in the counter-jet is 
the apparent radial velocity separation of the two knots. Us- 
ing the nomenclature of Micono et al. (1998), we identify the 
peak closer to the IR source as "c-jet a" which we will re- 
fer hereafter as the Counter-jet North, and "c-jet b+c+d" 
as the jet formed by three knots (unresolved in our images) 
which correspond to the southern knot and we will refer to 
it as the Counter-jet South. The brightest emission peak, 
Counter-jet North, appears at a distance of ~ 15" from 
HH46 with the emission centered at a heliocentric radial 
velocity of ~ +27kms~ 1 , suggesting that the counter-jet 
starts moving near the systemic velocity of the globule as 
it travels away from the source and gradually increases its 
radial velocity by « 28kms _1 within ~ 9" along the south- 
western emission lobe where the second peak, Counter-jet 
South, appears redshifted by ~ +55kms _ . We note also 
that this same pattern for acceleration of the gas is also 
seen in the proper motion data of Micono et al., where the 
proper motion knot shift is of order 0.082 arcsec yr _1 for 
knot c-jet a (Counter-jet North) and 0.138 arcsec yr _1 for 
knot c-jet d (farthest part of Counter-jet South). 

We have also detected some emission at the southern 
edge of the slit (row 1) which we identify as "knot 5" in 
the nomenclature of Eisloffel et al. (1994, see Fig. 2a) or 
"knot g" using the latest near-infrared mapping of this re- 
gion (Micono et al. 1998). We note that the position of this 
knot in the PV map suggests that the redshifted jet is flow- 
ing at essentially the same heliocentric radial velocity as the 
Counter-jet South knot. 

Although the counter-jet flow shows up very brightly in 
the H2 1-0 S(l) line emission map, it gets fainter with higher 
excitation energy of the H2 levels (see Fig. 3). Note how the 
counter-jet is just barely visible in the 2-1 S(2) emission 
map and completely disappears on the noise in the 3-2 S(3) 
emission map. The H2 2-1 S(l) contour emission map reveals 
that the counter-jet is still well defined, showing the same 
characteristic "double peak" shape as in the 1-0 S(l) map. 
Some fainter emission is still visible at the HH47A position 
but HH46 is now barely seen above the continuum emis- 
sion from the reflection nebula. The H2 1-0 S(0) emission 
map is very similar to the 1-0 S(l) map with HH46 stand- 
ing well out from the nebular continuum. This map also 



shows some residual OH lines that could not be perfectly 
subtracted from the final reduced object frame. Howe ver, 
this does not affect the flux measurements in Section 3.2 
since the OH lines are non- coincident with the H2 line. The 
emission maps obtained for the 3-2 S(l) and 4-3 S(3) lines 
(not shown) show substantially fainter emission from these 
objects implying a lack of high excitation H2 in the gas since 
these lines arise predominantly from hot spots in the flow 
where the temperature is high enough to populate signifi- 
cantly the upper energy levels of the H2 molecules. 

Given this, further analysis is thus restricted to the 
brightest four lines, 1-0 S(l), 2-1 S(l), 1-0 S(0) and 3-2 S(3). 
Table 3 shows the physical parameters of these H2 lines used 
in this work. 

In order to fully investigate the emission spectra from 
the counter-jet we coadded all the detector rows comprising 
the entire counter-jet flow (rows 10-20) and fitted simul- 
taneously two gaussians to the resulting spectra. Thus, the 
reference to Counter-jet North, in column 1 of Table 4, refers 
to the gaussian solution to the fit of the emission peak closer 
to the source whereas the reference to Counter-jet South 
refers to the gaussian solution to the fit of the emission peak 
far from the source. The reference to the Counter-jet Total 
was obtained by fitting a single gaussian profile to the spec- 
tra but allowing a broader FWHM to account for the line 
blending. The computed three-fold H2 line fluxes from the 
counter-jet are presented in Table 4. 

We have also searched for Brackett-7 emission from the 
counter-jet, but nothing was found. The 2.166 /im map only 
showed very weak continuum emission near the position of 
HH46 suggesting that the observed infrared emission is pro- 
duced mainly by dust scattering in the diffuse reflection neb- 
ulae associated with HH46. 



3.2 H2 excitation temperature 

The physical conditions for the excitation of H2 emission can 
be quantified by calculating the H2 level populations that are 
excited in the gas. Since H2 is a homonuclear molecule it may 
only be excited through electric quadrupolc radiation. In 
these conditions it is usual to assume the gas to be optically 
thin at these wavelengths, greatly simplifying the radiation 
transfer. Thus, the line intensity, I v ,j, for an H2 transition 
out of level (v, J), is directly proportional to the transition 
energy, hv V: j, the Einstein coefficient, A Vt j and the column 
density N Vt j of molecules excited in that particular level, 

Iv,j = hi> v ,jA Vl jN V: j/4ir 

where h is Planck's constant and v is the transition fre- 
quency. The values for A Vl j are those from Turner et al. 
(1977) and the values for v v% j are taken from Dabrowski 
(1984). For a Boltzmann distribution of the vibration- 
rotation levels (thermal distribution) at some local excita- 
tion temperature T ex , the population density, N V ,J, of level 
(v,J), is proportional to the statistical weight (level degen- 
eracy) g V: j and the Boltzmann factor. Thus, 



N v . 



j = g v ,j e 



where E v ,j is the transition energy, k is Boltzmann's 
constant and T cx is the excitation temperature. The exci- 
tation temperature can now be estimated from the data by 
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plotting N Vt j/g Vt j versus the upper level energy of the tran- 
sition E v j. The coefficient of the exponential fit to this col- 
umn density excitation diagram is proportional to T^" 1 . 

Figure 4 shows the H2 excitation diagrams calculated 
for the counter-jet emission using the data in Table 4. 
Since the counter-jet appears spatially extended and con- 
stituted by two emission peaks with different velocity com- 
ponents, we have plotted in Fig. 4 three separate data sets 
corresponding to the total flux emission from the whole 
counter-jet (Counter-jet Total), the northern peak (Counter- 
jet North) and southern peak (Counter-jet South) emission. 
The excitation temperature, T 2i i, indicated by the 1-0 S(l) 
to the 2-1 S(l) line ratio in the Counter-jet North position 
is 2100 K while in the Counter-jet Total position it is 1800 
K (see Table 4). A lower value of 1300 K is found for the 
Counter-jet South position. In Figure 4, the solid line repre- 
sents a least square fit through all the data points, providing 
an excitation temperature of T ox = 2600 K for the H2 gas at 
both the Counter-jet North and Total positions. The typical 
error of these estimates is given by the least squares fit and 
is of the order of 300 K. 

Also shown in Table 4 are the H2 line intensities and cor- 
responding column densities for the HH46 and the HH47A 
objects. Unfortunately, at these positions the low signal-to- 
noise obtained in the higher excitation transitions (v <; 3) 
hampers a more accurate determination of the excitation 
temperature. The excitation temperature for these objects 
are T 2 ,i=1200 K and T 2 ,i = 1800 K for the HH46 and HH47A 
objects, respectively. 

3.3 The [Fen] emission 

The forbidden emission from [Fen] at 1.257 /im (a 4 D 7 / 2 — 
a 6 D 9 / 2 ) and 1.644 /im (a 4 D 7 / 2 — a 4 F 9 / 2 ) was observed on 
the third night at the same slit position used for the H2 line 
observations to seek higher excitation regions within the flow 
and also to provide an estimate of the foreground extinction 
in the emitting gas. The position-velocity emission maps ob- 
tained for these two transitions are shown in Figure 5. Su- 
perposing both maps on a constant velocity grid makes all 
the bullet features coincide to within half pixel which gives 
an indirect way of checking the spatial shift correction per- 
formed on these maps. Since these [Fe n] PV maps show im- 
portant morphological and kinematical differences from the 
H2 maps presented previously, we will next examine them 
in detail. 

In the northeastern blueshifted outflow region, we see 
from Fig. 5, that this part of the flow appears composed of 
four bright [Fell] knots with V/fei ~ —250 to — 280 kms -1 , 
from which the northern most is clearly (spatially) associ- 
ated with the HH47A object (at the top edge of the array 
in the [Fen] 1.644 /urn map due to the spatial shift occurred 
when changing grating positions). This confirms the detec- 
tion made by Eisloffel et al. (1994) of a strong near-infrared 
[Fen] 1.644 /im emission peak at the position of the opti- 
cal HH47A object. We note that the radial velocity of the 
emitting [Fen] in the HH47A object (— 250 kms -1 ) is much 
greater than the corresponding H2 emission at — lOOkrns -1 . 
On the other end, and at the base of this blueshifted outflow, 
the HH46 object shows up brightly against a very weak con- 
tinuum. However, the most remarkable features in this Fig- 
ure, are the appearance of two emission [Fe 11] knots in the 



northeastern part of the HH 46/47 jet, in the region bridg- 
ing HH46 to HH47A, contrasting sharply with the morphol- 
ogy presented by the H 2 emission maps (Figs. 2-3). These 
two knots appear closer to HH46 and their locations in the 
outflow seem to have some parallel with the optical emis- 
sion. Several high-resolution optical images of this part of 
the outflow have been published that show detailed struc- 
ture across the whole region, from the base of the outflow at 
HH46, passing through the bright HH47A knot and ending 
in the HH47D object (Reipurth and Heatcote 1991, Harti- 
gan et al. 1993, Eisloffel and Mundt 1994, Heathcote et al. 
1996). The [Fen] maps presented here show the existence 
of two relatively bright knots in this part of the outflow 
which we can tentatively identify with the optical knot B8 
located at a distance of about 18" from HH46 (Eisloffel and 
Mundt 1994), and knot B0 (HH47B), located at a distance 
of about 32" from HH46 (see Figure 6). This is in contrast 
with the faint H2 outflow which appears, in this region, to 
fade away at approximately the position of the B8 [Fen] 
knot. Thus, it appears that the [Fen] emission is tracing 
a higher velocity outflowing jet component with no rela- 
tion to the lower velocity H2 jet. The whole northern [Fen] 
outflow appears blueshifted at heliocentric radial velocities 
between V He l ~ -280 kms -1 (HH46, B8 and HH47B) and 
V H ei ~ -250 knis" 1 (HH47A). 

The [Fen] 1.644 /im emission was first observed in this 
object by Eisloffel et al. (1994, see Fig. 2a). Their emission 
map, however, shows a slightly different spatial distribution 
of [Fen] 1.644 /im emission, especially regarding the HH47B 
object. While the blueshifted knots associated with HH46 
and B8 coincide in both maps, knot HH47B was not de- 
tected in their contour map. Two reasons might explain this 
fact. Either our map is slightly more sensitive or the evolu- 
tion of the gas emission over time may have determined the 
fading of the emission knot. Either way, we note that both 
emission knots are clearly visible in our [Fen] 1.644 /im and 
[Fen] 1.257 /jm maps and are both above (2<r) the noisy sky 
background. Overall, the spatial morphology of our [Fen] 
1.644 /xm PV map closely matches the [Fen] 1.644 /im im- 
age map of Eisloffel et al. although, the emission in our im- 
age is clearly more extended and brighter along the main 
blueshifted outflow axis. 

Turning to the southwestern region, there is a redshifted 
bright knot which is clearly associated with the emission 
from the counter-jet. As seen earlier, this counter-jet emis- 
sion is elongated and marginally resolved in to two H2 knots. 
However, we now see the [Fen] emission from the counter- 
jet as a single redshifted knot at Vuei = +100 km s -1 . More- 
over, when we overlay the H2 and [Fe n] emission (see Figures 
6 and 7), we conclude that this [Fen] knot is associated with 
the Counter-jet North peak and that no [Fen] emission is 
seen from the Counter-jet South peak. Notice also how this 
[Fen] knot seems to appear marginally closer to the out- 
flow source than to the bright Counter-jet North H2 knot 
(c.f Figure 7) while it appears to be travelling substantially 
faster than the corresponding H 2 knot. 

3. 3. 1 Reddening 

The extinction can be calculated using the relative flux ob- 
served from the [Fe n] transitions since these two lines arise 
from a common upper level of the Fe + ion (a 4 D 7 / 2 ) and the 
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line ratio is therefore independent of the number of excited 
ions on that particular level state. Thus, their intensity line 
ratio, /1.257//1.644, should depend only on the intrinsic prop- 
erties of the Fe + ion, provided that the lines are optically 
thin. This requirement should be true since the critical den- 
sity for the [Fen] 1.644 /itm transition is n cr it ~ 3 x 10 4 cm~ 3 
(Oliva, Moorhood and Danzinger 1990). The differential red- 
dening between the two [Fen] lines, Ej-h, can be computed 
from: 

/l.257//l.644 = (/l. 257 //l.644) X 

where (7i.257//i.644) is equal to 1.36; the theoretical 
line ratio given by the Einstein spontaneous probabilities 
of Nussbaumer & Storey (1988). The visual extinction is 
obtained using Eb-v ~ 3 x Ej-h (Draine 1989) and as- 
suming a normal reddening law: Ay = 3.1 x Eb-v- The 
derived [Fell] line flux and extinction estimates are shown 
in Table 2 for the observed emission peaks. This table shows 
that the [Fen] bullet associated with the Counter-jet North 
has a visual extinction of 7.05 ± 2.26 whereas an extinction 
value of 9.38 ± 1.49 is measured in the bullet associated with 
HH46, closer to the obscured YSO. Given the errors, these 
values indicate that the extinction is uniform throughout 
the region. 

Note also that we measure Ay ~ for the northern 
[Fell] bullets B8 and HH47B; given the uncertainties asso- 
ciated with the low flux measured for these bullets. This 
strongly suggests that the visual extinction is negligible in 
this (blueshifted) northwestern part of the outflow. However, 
as noted by Oliva and Origlia (1992), the [Fell] 1.644 /jm 
transition sits close (~ 40kms~ 1 ) to the relatively bright 
OH (5,3) -Ri(2) transition which hampers a more accurate 
determination of the reddening. 



4 DISCUSSION 

The intensity ratio of the w=l-0 S(l) line to the v=2-\ S(l) 
line has been widely used to discriminate between two com- 
mon types of physical mechanisms for the excitation of HH 
objects: thermal (shocks) and non-thermal (e.g. fluorescence 
by Lyman-a or UV continuum pumping, H2 reformation). 
HH objects are preferentially shock excited but fluorescent 
excitation by Lyman-a pumping may also be present (HH43; 
Bohm, Scott and Solf 1991, HH47A; Curiel et al. 1995, 
HH7; Fernandes and Brand 1995). The 1-0/2-1 S(l) ratio 
derived from the observations (Table 5) for the Counter-jet 
Total, Counter-jet North and HH47A objects are 15.7 ±2.3, 
10.1 ±3.5 and 16.4 ±2.5, respectively. These values are con- 
sistent with shock excitation (in which the 1-0/2-1 S(l) ratio 
<; 5) as fluorescent excitation requires the 1-0/2-1 S(l) ratio 
to be 2 in low density regions. In high density regions the 
fluorescence 1-0/2-1 S(l) ratio may be S> 2. However, this is 
unlikely to be the case here, since these high density models 
predict large column densities in the v=3 and v=4 levels 
and thus strong emission from lines such as the 3-2 S(3) and 
4-3 S(3) which could not be observed from the data. Under 
these conditions, the molecular hydrogen emission has its 
origin behind dense shocks and the H2 excitation diagrams 
calculated for the HH objects (Figure 4) are consistent with 
thermalised gas at a temperature of T^,i « 1300-2100 K for 



the counter-jet whereas HH46 and HH47A display T2,i ~ 
1200 and 1800 K, respectively. We note here that non-LTE, 
low density gas, may also produce line ratios similar to those 
observed here, as we will discuss below. 

The densities required to thermalise the H2 levels up 
to 20 000 K are of order 10 6 -10 7 cm~ 3 in a pure molecular 
gas (Burton et al. 1989). These high densities can be at- 
tained in regions close to the source from where the counter- 
jet emerges or in the dense relaxation medium behind the 
shock front where the gas cools via atomic, ionic and molec- 
ular line emission. Both pure hydrodynamic shocks (J-type) 
and magneto-hydrodynamic shocks (MHD C-type) can pro- 
duce significant emission from H2, although only fast (dis- 
sociative) J-type shocks can produce observable [Fe 11] emis- 
sion since C-shocks attain a lower post-shock temperature 
2000 K) than J-shocks (> 2000 K). 

Comparing the 1-0/2-1 S(l) line ratios with the shock 
models of Smith (1995) and Eisloffel et al. (1996) yields a 
useful constraint on the various critical shock parameters 
such as the pre-shock density and the type and geometry of 
the shock. First, note from Table 5 that the observed range 
for the line ratios of 10-17 (except for the noisier Counter-jet 
South and HH46 objects) is consistent with a slow planar J- 
shock with speed 8-9 km s -1 . However, these line ratios can 
be also reproduced by fast planar C-shocks with speeds of 
35-40 km s" 1 . The modelled 1-0 S(l)/3-2 S(3) line ratio of 
447.7 for a 35 kms^ 1 C-shock is nonetheless much greater 
than the observed value of only 45-50 and thus planar C- 
shocks can be ruled out. On the other hand, the slow J- 
shocks are much weaker than those required to excite the 
observed [Fell] lines as we will see below. On these grounds, 
we conclude that neither planar type of shock can explain 
the observations, except, perhaps, for the Counter-jet South 
and HH46 objects from which no conclusions can be drawn 
due to the large errors in the observed line ratios. 

Constant temperature models were also compared with 
the observations. The single (IT) and mixed temperature 
(2T) slab of gas models for both LTE (n=10 s cm~ 3 ) and 
non-LTE conditions (n=10 4 cm -3 ) yield specific line ratios 
which are given in Table 5 (Smith, private communication). 
The three chosen IT models in Table 5 fail to account for 
both line ratios regardless of the LTE conditions. Moreover, 
the excitation temperature derived from the 1-0/2-1 S(l) 
line ratio ranges 1800-2100 K, whereas the least square fit 
to all the data points implies a higher excitation temperature 
of 2600 K. Given this, we may expect better results from the 
dual temperature models. In fact, an inspection of Table 5, 
clearly indicates a closer match to the data. However, there 
are strong physical constraints on these approach since we 
would need a heating mechanism capable of maintaining the 
gas at fixed temperatures of order of a few thousand Kelvin 
for a timescale of years. At present, there is no valid expla- 
nation for such mechanism, and thus we also rule out this 
possibility. 

In agreement with the observed morphology for these 
type of objects, we thus turn our attention to the exci- 
tation of both H2 and [Fell] emission in bow-shock mod- 
els. Strong [Fell] emission is often interpreted as indicative 
of shock excitation (e.g Moorwood & Oliva 1988) and can 
be expected in, for example, supernova remnants (McKee, 
Chernoff &: Hollenbach 1984). In bow-shocks, however, the 
high excitation [Fell] emission derives from the bow-head, 
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where the shock is strongest (see e.g. HH1; Davis, Eisloffel 
& Ray 1994). Further downstream, in the lower excitation 
oblique bow wings, emission from H2 can be excited in the 
softer (non-ionising) shocks. The presence of [Fell] 1.644 /im 
in the flow is thus considered a diagnostic of fast dissociative 
shocks at the tip of the bullet that heads a post-shocked and 
accelerated medium observed as a trailing H2 wake (see e.g. 
Tedds, Brand and Burton 1998). 

The bow shock emission structure can be modelled by 
a surface of revolution which can be simply written as 2 oc 
r s in cylindrical coordinates, where s is the free parameter 
that adjusts the bow aperture (see Fig. 6 of Eisloffel et al. 
1996). This morphology combines in a single structure, two 
distinct regions, depending on the local physical conditions, 
especially the ion fractions and the strength of the magnetic 
field. Thus, we need to distinguish between J-type or C-type 
bow shocks. 

The former, produces emission which is everywhere ex- 
cited by pure hydrodynamic type shocks, from the apex of 
the bow, where strong J-shocks destroy the H2 molecules 
and ionises Fe, to the extended H2 emitting wings down- 
stream along the surface. The predictions from these bows 
depend critically on the shape parameter s. The best param- 
eter fit models are shown in Table 5. In general, both LTE 
and non-LTE models come close to adjust the observed line 
ratios. The best fit bow shape parameter in these models 
range from 1.2 to 2.0 (paraboloidal) which is consistent with 
strong [Fell] emission over much of the outflow. Narrow line 
profiles are expected from these type of shocks and can be 
used to discriminate between shock types, since the larger 
shock front of C-type shocks will produce wider line profiles. 
This, however, requires high-resolution spectroscopy data 
and therefore with the present set of data we cannot rule 
out the J-shock bow model for the Counter-jet or HH47A. 

On the other hand, the C-type bow shock model dif- 
fers from the above picture only on the bow wings, where 
the softer C-type shocks produces emission from H2 and CO. 
The Orion bullets are the best studied example of this struc- 
ture, where a high- velocity leading [Fell] bullet is plough- 
ing through the medium leaving behind a wake of lower- 
excitation H2 emitting gas generated at the warmer and 
oblique C-type shocks at the bow flanks (Tedds, Brand and 
Burton 1999). The C-type bow-shock model fits shown in 
Table 5 can also be made to agree with the present line ra- 
tios, but for a larger s bow shape parameter (1.5 < s < 2.0). 
In this case, the bows are closer to paraboloids and less 
pointed than the J-type bows. 

The [Fen] 1.644 /im surface brightness predicted for a 
fast (v > 30kms -1 ) J-shock can be large due to the large 
column of warm gas generated by UV absorption and H2 
reformation in the post-shock layer (Smith 1994). For a 
100 kms -1 shock with a pre-shock density of 10 4 cm -3 , a 
surface brightness of 1.5 x 10 -4 ergs -1 cm -2 sr -1 for the 
[Fen] 1.644 /im line can be achieved. This value is consis- 
tent within a factor of 3 with our measured value at the 
Counter-jet North peak of 4.6 x 10 _4 ergs -1 cm -2 sr -1 . Note 
that this is as strong as the observed H2 1-0 S(l) surface 
brightness of 4.4 x 10 -4 ergs -1 cm -2 sr -1 . Such high values, 
where the [Fell] 1.644 fim /H2 1-0S(1) line ratio are of equal 
strength, have already been observed in other HH objects 
by Stapelfeldt et al. (1991) and can be well modelled with 
fast J-shock models (Smith 1994). 



A non-dissociative J-shock (v ^ 22 kms -1 ), on then 
other hand, produces smaller columns of warmer gas. We 
may estimate an upper limit for the [Fen] 1.644 /im line 
strength from the model employed by Smith (1994) under 
the assumptions that the electron fraction is ~ 10 -3 and the 
abundance of Fe + is ~ 10 -6 . For the fastest non-dissociative 
shock (with v = 22 kms -1 ), the line strength predicted is 
< 4.2 x 10 -5 ergs -1 cm -2 sr -1 , a value which is about 10 
times smaller than observed. For a 15 kms -1 shock speed, 
the [Fell] 1.644 /im intensity drops down by a factor of 2 
and therefore we rule out the slow J-shock model for the 
generation of excited [Fell]. 

Downes (1997) has modelled the H 2 and [SB] A6731 
emission from the HH 46/47 counter-jet using an hydrody- 
namical code to predict intensity maps and line profiles. 
The simulated emission showed the appearance of numer- 
ous mini-bowshocks throughout the outflowing jet result- 
ing from Rayleigh- Taylor instabilities. Moreover, the synthe- 
sised [SIIJA6731 and H2 1-0 S(l) line profiles derived from 
their emission maps of the counter-jet, yield clear asymme- 
tries in the profiles with the near-infrared line peaking at 3- 
4kms -1 , while the optical line peaks at about 150 kms -1 . 
Since the H2 emission appears only at the bow flanks, it 
displays a low velocity component; in our data +27 and 
+55 kms -1 , for the Counter-jet North and South objects, 
respectively. In contrast, the [Fen] knot coincident with the 
Counter-jet North object, is seen peaking at a higher ve- 
locity of + 100 kms -1 . This agrees well with the bow shock 
scenario, where the high excitation lines of [Fell], [SII] and 
Ha, are produced in the strong and high speed shocks near 
the bow apex, yielding higher normal gas velocities of ~ 100- 
150 kms" 1 . 

Looking at the relative position of the [Fell] and H2 
emission in HH47A (see Figs. 6 and 7), we see the [Fell] peak 
ahead about 5" from the H2 peak. The spatial calibration 
is accurate to less than half pixel (~ 1") and thus we con- 
sider the shift to be genuine. This translates to a separation 
distance of ~ 3 x 10 16 cm from the [Fell] emitting section 
to the H2 emitting section. This separation is characteris- 
tic in bow-shocks viewed at orientation angles > 30 degrees 
to the line of sight (Smith 1991) and thus we suggest that 
HH47A can be well modelled by a bow shock working sur- 
face, in a similar fashion to HH99B (Davis, Smith & Eisloffel 
1999). Note, however, the large radial velocities achieved for 
the [Fen] (-250 kms" 1 ), Ha (-130 kms" 1 ) and even H 2 
(— 100 kms -1 ), whose molecules should be all dissociated at 
these high (v J> 50 kms -1 ) speeds. To produce H2 emission 
at these extreme velocities, it is necessary that the pre-shock 
medium had already been set in motion by some previous 
event (e.g. through the passage of an earlier shock wave 
which may be accounted by the faint HH47D object), so 
that the relative velocity between the two flows stays lower 
than the dissociative limit of the H2 molecule: 25 kms -1 , 
for a J-shock or 50 kms -1 , for a C-shock (Smith 1991). Un- 
fortunately in this work, we have not reached the required 
sensitivity in the 3-2 S(3) line in HH47A to help discriminate 
between these two types of shocks spectroscopically. Thus, 
on geometrical grounds, we expect HH47A to be a classical 
bow shaped shock, with a dissociative J-shock section at the 
bow apex (which produces the [Fell] and the optical emis- 
sion), plus a longer C-shock wing section which produces the 
lower excitation H2 emission. With the measured tangential 
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velocities in this object of the order ~ 180 km s _1 (both opti- 
cal and near- infrared) , a bow shock velocity of w 220 kms -1 
can be inferred. 

The emission from HH46 is partially blended with the 
reflection nebula which causes the Eh emission to be scat- 
tered and thus appear weaker in contrast to the brighter 
[Fell] emission. The 2-1 S(l) line emission from HH46 is 
somewhat below our expectations since the 1-0 S(l) re- 
mains strong. In fact, the 2-1/1-0 S(l) line ratio is only 
0.014 compared with, for example, 0.061 in HH47A. We 
suggest that we may be observing only the emission tail 
of HH46, which produces weaker shocks, and thus explain- 
ing the weaker H2 2-1 S(l) emission observed and also the 
— 70kms~ heliocentric velocity shift observed for this ob- 
ject. Moreover, the [Fen] emission peaks at — 280kms~ 1 , 
implying very fast bow shock speeds near the source. In 
fact, a shock which produces [Fell] efficiently, will likely de- 
stroy the H2. For H2 to survive the extreme shock velocities 
observed in HH46, we therefore require a mechanism that 
can soften the shock. Commonly cited in the literature are 
mechanisms such as the magnetic precursor model, the en- 
trainment of material from the wind/ambient material inter- 
face or the high Alfven Mach number shock absorbers (bow 
shocks with a high magnetic field). Distinguishing between 
these physical processes is beyond the scope of this paper 
with the available set of data. Obtaining higher resolution 
data will certainly help to establish the importance of these 
mechanisms. 



5 CONCLUSIONS 

The physical conditions along the redshifted counter-jet 
and the blueshifted HH46 to 47A outflow were investigated 
through observation of velocity-position diagrams and spec- 
tra in the J ([Fell] 1.644 /im and [Fell] 1.257 (im) and K 
(several H2 transitions) infrared bands. 

The observed H2 line emission confirms the existence of 
a bright and extended redshifted counter-jet outflow south- 
west of HH46. Moreover, there seems to be compelling ev- 
idence for acceleration of the counter-jet from 27kms~ 1 to 
55kms _1 along the southwestern emission lobe as it travels 
away from the source. 

The observed [Fen] emission suggests a near uniform 
extinction in the HH46 and the counter-jet region of w 7-9 
visual magnitudes, whereas the extinction calculated for the 
blueshifted knot B8 and HH47B objects is negligible. 

The H2 excitation temperatures derived from the 1-0/2- 
1 S(l) ratio are in the range 1800-2100 K for all the objects, 
with the exception of the HH46 object for which the excita- 
tion temperature is 1200 K. The unusual low 2-1 S(l) line 
flux in this object was probably caused by the limited region 
covered by our slit in this object which may have sampled 
only the wing sections of the bow-shock structure. The in- 
clusion of the 3-2 S(3) in the H2 excitation diagrams for the 
counter-jet shock yields an excitation temperature for the 
gas of T ox ?a 2600K. 

The non-detection of higher energy transitions of H2 in 
the K band, precludes several non-thermal excitation mecha- 
nisms such as fluorescence or H2 reformation. The observed 
2-1/1-0 S(l) and 3-2/1-0 S(l) line ratios were found con- 
sistent with models of J and C-type bow shocks, whilst 



ruling out planar shocks. In the counter-jet, the shape pa- 
rameter for C-type bows can be constrained to the range 
1.5 < s < 2.0 whereas J-type bows need to be made ex- 
tremely wide, 1.2 < s < 1.5, to adjust well our data. 
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Table 1. 


Parameter list of the observed 


emission 


lines 








An 


Total Exj 


3. Spatial 


Resolution 


Species 


label 


t M m ) 


(min)r 


Shiftl] 


Av (kms" 1 )! 


H 2 


1-0 S(l) 


2.1218 


3 


0.0 


66.8 


H 2 


1-0 S(0) 


2.2235 


1 


1.0 


61.5 


H 2 


2-1 S(l) 


2.2477 


8 


1.5 


60.3 


H 2 


2-1 S(2) 


2.154 


5 


0.5 


65.1 


i_i 

ti2 


2-1 S(3) 


2.065 


Q 
O 


n e: 
-U.o 


/ U.U 


H 2 


3-2 S(l) 


2.3846 


8 


2.5 


53.6 


H 2 


3-2 S(3) 


2.2014 


10 


0.5 


62.6 


H 2 


4-3 S(3) 


2.3445 


8 


2.5 


55.5 


H 


Brackett-7 


2.166 


5 


0.5 


64.5 


[Fen] 


a 6 Dg/ 2 - a A D 7/2 


1.257 


2 


3.5 


47.7 


[Fen] 


a 4 F 9/2 - a 4 D 7/2 


1.644 


1 


-2.5 


97.9 


a This is 


the total exposure 


time of the final reduced frame. 





b Spatial shift or displacement (in pixels) between frames due to the use of different grating settings which were calculated 
using observations of the standard star Eta Canis-Majoris at different wavelength settings. 
c Velocity resolution per pixel calculated using the wavelength resolution: Av/v = A A/ A 

Table 2. Line flux and extinction estimate towards HH46, the counter-jet and the [Fell] bullets 



Parameter 


Counter-jet North 


HH46 


knot B8 


HH47B 


[Fen] 1.257 /an 
[Fen] 1.644 /an 


3.95 (0.87) 
5.84 (0.22) 


4.85 (0.67) 
9.03 (0.45) 


4.79 (0.77) 
3.10 (0.30) 


3.98 (0.87) 
2.44 (0.34) 


1.257/1.644 

A2.1218 
E s -v 
Av 


0.68 (0.15) 
0.76 (0.24) 
0.78 (0.25) 
2.27 (0.73) 
7.05 (2.26) 


0.54 (0.08) 
1.01 (0.16) 
1.04 (0.16) 
3.03 (0.48) 
9.38 (1.49) 


1.54 (0.29) 
-0.14 (0.20) 
-0.14 (0.21) 
-0.41 (0.61) 
-1.28 (1.90) 


1.63 (0.42) 
-0.20 (0.28) 
-0.20 (0.29) 
-0.59 (0.84) 
-1.83 (2.61) 



Note: The line flux are in units of 10 



Table 3. Physical parameters for the brightest H 2 lines 

Transition 
A (/an) 
Energy (Kelvin) 
Statistical weight 
Spontaneous decay rate 
Extinction (Ax) 



1-0 S(l) 


1-0 S(0) 


2-1 S(l) 


3-2 S(3) 


2.1218 


2.2235 


2.2477 


2.2014 


6956 


6471 


12550 


19086 


21 


5 


21 


33 


3.47 


2.53 


4.98 


5.63 


1.04 


0.96 


0.94 


0.98 



a The statistical weight is (2J+1)(2/+1); where J is the rotational quantum number and I is the spin quantum number 
b Einstein coefficient in units of 10 -7 s~ 1 



© 0000 RAS, MNRAS 000, 000-000 



2 



Table 4. Observed H2 line intensities and derived column densities 



Object 


Transition 
label 


Observed 
line intensity^ 


De-reddened 
line intensity^ 


(cm-2)| 


qpex 
2 1 

(Kelvin)^ 


(Kelvin)[] 


HH46 


1-0 S(l) 

1- S(0) 

2- 1 S(l) 

3- 2 S(3) 


280 (32) 
33.5 (6.6) 
4.28 (2.25) 
< 21.8 


57.4 (6.6) 
6.5 (1.3) 
0.82 (0.43) 
< 4.3 


105.8 (12.1) 
72.3 (14.3) 
1.12 (0.59) 
< 3.2 


1200 




HH47A 


1-0 S(l) 

1- S(0) 

2- 1 S(l) 

3- 2 S(3) 


376 (29.9) 
29.9 (5.5) 
24.5 (3.1) 
< 25.1 


77.1 (6.1) 
5.8 (1.1) 
4.7 (0.6) 
< 4.9 


142.1 (11.3) 
64.6 (11.9) 
6.4 (0.8) 
< 3.7 


1800 




Counter-jet 
North 


1-0 S(l) 

1- S(0) 

2- 1 S(l) 

3- 2 S(3) 


212.9 (58.5) 

37.4 (5.8) 

22.5 (4.8) 
5.0 (2.3) 


43.7 (12.0) 
7.3 (1.1) 
4.3 (0.9) 

0.98 (0.45) 


80.4 (22.1) 
80.7 (12.5) 

5.9 (1.3) 
0.74 (0.34) 


2100 


2600 


Counter-jet 
South 


1-0 S(l) 

1- S(0) 

2- 1 S(l) 

3- 2 S(3) 


138.3 (58.9) 
17.5 (4.45) 
2.91 (2.6) 
< 18.0 


28.4 (12.1) 
3.40 (0.86) 
0.56 (0.50) 
< 3.5 


52.3 (22.3) 
37.8 (9.6) 
0.76 (0.68) 
< 2.6 


1300 




Counter-jet 
Total 


1-0 S(l) 

1- S(0) 

2- 1 S(l) 

3- 2 S(3) 


364.1 (23.5) 
54.6 (3.0) 
24.8 (3.2) 
7.9 (1.1) 


74.7 (4.8) 
10.6 (0.6) 
4.8 (0.6) 
1.6 (0.2) 


137.6 (26.6) 
117.9 (19.4) 
6.5 (2.5) 
1.2 (0.5) 


1800 


2600 


a Units of 10" 


~ 6 ergs -1 cm 


- 2 sr-! 











b Extinction corrected line intensities in units of 10 -5 ergs -1 cm -2 sr _1 according to the values in Table 2 

c Integrated column density divided by the number of degenerate levels (statistical weight) in units of 10 13 cm -2 

d Excitation temperature calculated using the 1-0 S(l) to 2-1 S(l) line ratio 

e Excitation temperature calculated using a least square fit to all the data points 
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Table 5. Line ratios compared with planar J and C-type shock models from Smith (1995). 



Object or JVIodel 


h-oS(i)/h-iS(i) 


Il-0S(l)/l3-2S(3) 


HH46 


70.0 (37.6) 


> 13.4 


HH47A 


16.4 (2.5) 


> 15.6 


Counter-jet North 


10.1 (3.5) 


44.5 (23.8) 


Counter-jet South 


50.8 (50.3) 


> 8.0 


Counter-jet Total 


15.7 (2.3) 


48.1 (7.4) 








Model J8 


17.9 


340 


Model J9 


10.0 


88 


Model J15 


4.2 


12.6 


Model J22 


3.51 


< 10 


Model C35 


20.1 


447.7 


Model C40 


10.9 


122 


Model C50 


4.75 


20 


Constant single (IT) and dual (2T) temperature slabs 






Model IT 1 


10 


112 


Model IT 


15 


277 


Model IT 


6.5 


45 


Model 2T 2 


10 


50 


Model 2T 


15 


100 


Model 2T 


9.6 


45 


Bow-shocks (J-type) 3 






Model LTE 


10 


56 


Model LTE 


15 


108 


Model LTE 


8.8 


45 


Model non-LTE, 2 X 10 4 cm" 3 


10 


48 


Model non-LTE, 2 X 10 4 cm~ 3 


15 


91 


Model non-LTE, 2 X 10 4 cm~ 3 


9.6 


45 


Bow-shocks (C-type) 4 






Model 10 cm ° 


10 


56 


Model 10 6 cm~ 3 


15 


111 


Model 10 6 cm- 3 


8.7 


45 


Model 10 4 cm- 3 


10 


50 


Model 10 4 cm~ 3 


15 


100 


Model 10 4 cm~ 3 


9.3 


45 


Model 10 8 cm~ 3 


10 


40 


Model 10 8 cm~ 3 


15 


67 


Model 10 8 cm~ 3 


11 


45 



1 LTE standard conditions, n=10 8 cm 3 (idcntical results are found for non-LTE conditions, n=10 4 cm 3 ). 

2 Mixed 1600K (95%) plus 3500K (5%) gas. 

3 Bow shape parameter: 1.2 < s < 1.5 

4 1.5 < s < 2.0 
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